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Abstract. 

In this work, we inspect the reliability of controlling and quelling an epidemic 
disease mimicked by a Susccptiblc-Infcctcd-Susccptiblc (SIS) model defined on a 
complex network by means of current and implcmcntable quarantine and isolation 
policies. Specifically, we consider that each individual in the network is originally 
linked to two types of individuals: members of the same household and acquaintances. 
The topology of this network evolves taking into account a probability q that aims 
at representing the quarantine or isolation process in which the connection with 
acquaintances is disrupted according to standard policies of control of epidemics. 
Within current policies of self-isolation and standard infection rates, our results show 
that the propagation is either only controllable for hypothetical rates of compliance or 
uncontrollable at all. 
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1. Introduction 

Alongside with wars and natural disasters, plagues and epidemic (pandemic) 
diseases lie at the top of death toll lists in Human history. By directly interfering 
with the risk of death, such events or even their likelihood are the source of acute 
distress among populations, which has been well-documented by historians since ancient 
times p]. Concerning plagues and epidemic diseases, isolation of infected people has 
been applied, at least, since Old Testament period as an instrument for controlling and 
quelling the spread of both viruses and contamination agents in such situations [2]. If 
for centuries it was possible to estrange people from their relatives and whereabouts by 
means of a simple decree, social conquests based on the "Declaration of the rights of 
man and of the citizen", which led to "Universal declaration of Human rights", have 
urged the discussion over the morality and legitimacy of compulsory quarantine and 
isolation (Q&I) practices [31 [U El [TJ. Indeed, countries like Brazil and Japan have 
repealed previous laws for compulsory isolation, abortion and sterilisation of patients 
suffering from leprosy as well as decided to pay compensations and allowances to people 
subject to these practices 0, EJ HQ]. Accordingly, current Q&I practices strongly rely 
on the consciousness of the infected individual and her attitude towards the rest of the 
society by imposing a self-isolation spell according to a medical recommendation (see [7] 
and references therein). Concomitantly, recent polls have shown that the willingness 
to comply with a self-isolation period depends on the social condition and literary 
habilitation of the individual very much [TT] . 

Quantitatively, the problem of reasoning over the spreading of an infectious and 
epidemic disease is generally based on standard models such as the Susceptible-Infected 
(SI), the Susceptible- Infected-Susceptible (SIS) and the Susceptible- Infected- Recovered 
(SIR) models [12], [13] and their variants. Despite the simplicity of these models, they 
have been successfully applied to a variety of cases [HI HHJ [161 El EE! EES EQ~1 [21] in 
the assessment of the propagation of an epidemic disease. In the form of differential 
equations, their dynamics is well-known in the literature of the theme where they 
are usually called general solutions. When geometrical conditions are taken into 
account, namely the structure of the network of relationships between people, the 
critical behaviour changes [22], [23] . Within this context, analytical and numerical results 
concerning the existence or not of a non-zero transition for quenched scale-free networks 
has been in the spotlight [21 ESI ESI E3 EE1 EHJ [30] . 

In our case, we have opted to base our model on the SIS model, which can 
represent situations wherein the virus responsible for the infection is able to mutate 
fast. Explicitly, although someone is recovered from the infection, she is at risk of being 
infected once again due to the newest mutation. Concerning the topology in which the 
phenomenon endures, we have assumed a dynamical complex network with N nodes, 
each one linked to k other nodes in accordance with a distribution P(k). Each node 
of the network represents an individual either in a Susceptible, S, or in an Infected, 
/, state. As usual, we start the dynamics with a single infected individual, randomly 



Probing into the effectiveness of self-isolation policies in epidemic control 



3 



chosen among the N nodes. At first, all the remaining individuals are Susceptible. In 
addition, we have considered that a density of the connections, d, in the network are 
fixed and thus they cannot be removed during the Q&I process. Two people connected 
by a fixed link are defined as members of the same household, whereas people linked by 
non-fixed links are dubbed acquaintances. 

Recently, some works have studied the effect of isolation in the physics and computer 
science literature [3H E21 E31 EU EH [361 EZ] • It was therein argued that such policies can 
be effective in controlling the spread of an epidemic disease. However, the conditions of 
these models rely on the total isolation of the infected patient, which is closer to academic 
proposals that also inspired several literary contemporary masterpieces [381 EH] rather 
than an actual implementable policy following international rules or recommendations. 
Alternatively, total isolation scenarios have also been studied in the context of small 
islands. In this case, it is argued that these policies might be successful, but in quite 
stern conditions demanding a complex and expensive logistics [4*0| WT\ . On the other 
hand, the study of the impact of the household in the propagation of SIR epidemics as 
well as adequate vaccination policies have been discussed in a biological and medical 
framework [121 E31 HU HH1 SHI EZ]- However, for the sake of analytical treatment [46J, 
they often neglect the structure of social networks and assume a random contact 
approach between individuals. Complementary, surveys stemming from the analysis 
of data on fatalities due to the "Spanish flu" pandemic in cities of the United States 
of America and its relation to public health and non-pharmaceutical interventions have 
been presented [HI H9] . In spite of the fact that both results provide important insights 
into the relevancy of coordinated interventions and respective set about, (American) 
society has dramatically changed in almost 100 years and thus its epidemic response 
to such interventions. Moreover, some of the measures adopted during that pandemic 
are now liable of being judged unethical or breaching some fundamental law, as we 
have previously mentioned to. Last but not least, it was impossible in these studies to 
separate out the impact of each measure. 

Thence, with this work we intend providing an answer to the fundamental question 
about to what extend Q&I policies abiding by World Health Organization (WHO) 
directives, which are prone to self-conscious isolation, are real hinderers of the spreading 
of an epidemic disease, contributing in this way to its control and suppression. 
Quantitatively the task of answering to this question is presented in the form of looking 
for an infection rate (or epidemic) threshold, A*, for a given rate of effectiveness of the 
isolation, above which the disease persists or it is quelled otherwise. 

2. Model and Numerical Results 

At each time step, t, the following automata rules control our model: 

• We visit every node in the network; 

• For each Infected node, we look for all Susceptible neighbours of hers; 
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Figure 1. Time evolution of the density of Infected individuals for d = 0.2 and q = 0.3 
(left side) and q = 0.8 (right side), and typical values of the infection probability A for 
random networks with N = 10 4 nodes and (k) = 5. Notice that the time needed for 
the system to reach the steady states increases for increasing values of the self-isolation 
probability q of the Q&I process. Results are averaged over 10 3 realisations. 



• For each I — S link connecting members of the same household, the Susceptible 
individual becomes Infected with probability A; 

• On the other hand, we remove each I — S link connecting acquaintance individuals 
with probability q (quarantine). Then, the Infected individual is isolated and the 
Susceptible one is reconnected to another randomly chosen Susceptible node, who 
is not a current contact of hers. For each surviving I — S link, the Susceptible 
individual becomes Infected with probability A; 

• After the verification of all neighbours of an Infected individual, she returns to the 
Susceptible state with probability a. 

Let us elaborate upon the parameters and the rules we have established. We start 
by discussing the role of q the main purpose of which is to reproduce the probability 
that a patient complies with the medical recommendation to stay home. Ideally, people 
would strictly abide by the physician's counsel, but as recent polls have showed and 
given statistical significance, the rates of compliance are different from 100% and can 
go as low as 50%. Two main factors for this change have been identified JTTJ [50]; i) 
More often than not people do not understand or distort what they are expected to 
do due to literacy skills [50]; ii) As occurs in life in general, people (patients included) 
perform a risk assessment [H] before making a decision on stopping contact with a 
given individual. Accordingly, infected agents decide to disrupt the contacts they have 
depending on the time they were asked to stay home, the type, namely the importance 
(either circumstantial or not) of the relationship between the infected person and 
the susceptible individual, and the number of people that can be involved when the 
attendance of a given event is considered. Furthermore, this risk assessment is weighed 
by the hazard that the eventual non-compliance carries. Effectively, the perception of 
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risk is quite affected as the time of the self-isolation goes by as well as the reckoned 
performance of pharmaceutical interventions. As a matter of fact, glaring differences in 
the polls results concerning self-isolation before and after the introduction of a vaccine 
in the H1N1 pandemics were verified. A similar behaviour has been verified in polls on 
high risk unprotected sexual relations as well [51] . 

Moving ahead, we now underpin the introduction of the rewiring process. As more 
than well known, social relationships come to pass because people need one another. 
Besides friendship, familiar and sentimental ties, people relate with other people because 
of their functions, skills and responsibilities. Previous studies have considered that 
after a link between susceptible and infected individuals having been destroyed, no new 
relation is established. In a sociological perspective, this option acts as though the needs 
of the susceptible individual can be solved within her current circle. Our option goes 
along another path we consider more realistic and that is close to processes of rewiring 
studied in other works [31], [52j [53] . In spite of the fact that we could have considered 
a probability / that a susceptible "stranded" agent relinks, we have opted to mitigate 
the number of parameters of the problem and keep it constant and equal to 

As the majority of people live together, in very few (and negligible) cases, it is 
possible to someone to really isolate herself from the people with whom she shares 
an address. In this way, we justify the introduction of parameter d. Still, we would 
like to shed light on some aspects of this parameter. First, its introduction does 
not transform our network into a superposition of two networks (householders and 
acquaintances) |§] The network remains unique, but there is a fraction d of the links 
that cannot be changed because they have been quenched. Householders are the people 
who share the same residence whether they are relatives or not (e.g., student's flats) 
and acquaintances are people with whom someone maintains contact with in a close and 
regular basis representing at its best a fraction of the so called Dunbar's number, which 
is an anthropological/cognitive measure that quantifies the number of people with whom 
each person is able to keep on a relationship, i.e., preserving some information about his 
life. It is also important emphasise that in establishing the network we have borne in 
mind important details, e.g., if two people are both members of the same household of a 
third individual, then a quenched connection between them is immediately established 
because they are obviously members of the same household. 

As regards other computational details, we have applied a synchronous update, 
i.e., after having visited all the nodes on the network, we update all the individuals 
states simultaneously [5H [55] . In other words, a new time step is considered only when 
the state of all the individuals has been scrutinised. Furthermore, we have considered 
in this work a fixed / — > S transition probability a = 0.2 and we have simulated the 
model on a Erdos-Renyi network with N = 10 4 nodes and (k) = 5. The results herein 

f We can assume that neither her householders nor acquaintances manage to act as a stand-in. 
§ From a coarse grained point of view, we could use this approach to define a network households but 
this is not a superposition of networks as well since in this case each node is now a household and not 
an individual. 
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Figure 2. Left side: Normalised stationary density of Infected individuals I s t as a 
function of A for d = 0.2 and typical values of q. We can observe a transition from 
I s t = to I s t > at different threshold points X(q). Right side: I s t as a function of 
the self-isolation probability q of the Q&I process, for d = 0.2 and typical values of the 
infection probability A. We also observe a transition for each value of A, at different 
points. In all simulations, N = 10 4 , (k) = 5 and averages are over 10 3 realisations. 

presented have been obtained averaging over 10 3 independent simulations, for each set 
values presented. Please observe we have only considered the epidemic spreading in the 
largest connected component of the network, i.e., the giant cluster [32] . 

In Fig. [U we exhibit the time evolution of the density of Infected individuals I(t) 
for 20% of fixed links in the network (i.e., for d = 0.2), two different values of the 
self-isolation probability, namely q = 0.3 and q = 0.8, and A scanning the domain of 
infection rates. As it can be observed, there are different values of A above which we have 
I st > 0. Here, the notation I st = lit — > oo) stands for the stationary density of Infected 
individuals jjj] This defines the usual phase transition of epidemic models: for A < A*(g), 
we have a disease-free phase with all the individuals presenting a S state, whereas for 
A > X*(q) we have an epidemic phase, i.e., the disease spreads out and a finite fraction 
of the population is constantly infected. Pay heed to the fact that, for increasing values 
of q, the time needed for the system to reach a stationary state increases. This kind 
of transition is best observed in Fig. [2] (left panel), wherein we exhibit the stationary 
density of Infected individuals I st as a function of A for d = 0.2 and typical values 
of the self-isolation probability q of the Q&I process. In this case, we can see that 
for a fixed value of q the transition occurs at different values of A, which defines the 
transition points A*. Analogously, the stationary density of Infected individuals I s t can 
be represented as a function of the self-isolation probability q of the Q&I process for 
different values of A (see Fig. [21 right panel). Once more, we perceive a transition at 
different values of q(X). Notice that when we increase the self-isolation probability, the 

|| The quantity I s t was obtained from time averages of 1(f) taken in the steady state, and in addition 
we also considered averages over 10 3 network realisations. 
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Figure 3. Normalized stationary density of Infected individuals I st as a function of A 
for d = 0.4 and typical values of q. We can also observe a phase transition from I st = 
to I st > at different threshold points X(q). In all simulations, N — 10 4 , (k) = 5 and 
averages are over 10 3 realizations. 



q a c x 2 

0.3 0.62 ±0.04 0.91 ±0.08 1.13 x 10" 4 

0.5 0.34 ±0.01 0.93 ±0.05 2.46 x 10~ 4 

0.8 0.11 ±0.01 0.73 ±0.09 6.28 x 10~ 4 



Table 1. Values of the parameters a and c used in Fig. |U 



disease only disseminates through the network if we increase the infection probability 
A. In other words, the final size of the epidemic may be reduced, as was observed in 
the 2009 H1N1 pandemic [H] and as also discussed in some works [3U |32} [331 EH 152] . 
Nonetheless, in what follows, we will conclude its outcome is limited when we take into 
account a more realistic model with a density of fixed links. 

We have simulated the model for other values of the density of fixed links d, and we 
have also observed the above-discussed transition (see Fig. [3]). Alternatively, we take 
into account the transition values A* as functions of d for fixed values of q and we have 
noted that the data obtained from the simulations are well described by a stretched 
exponential function, 

X*(d) ~exp[-(d/a) c ], (1) 

for different values of the parameters a and c (see Fig. H]). In particular, for small values 
of q like q = 0.3, the dependency of A* with d is almost purely exponential, since we 
have obtained an exponent c = 0.91 ± 0.08. The values of a and c are shown in Tab. 1. 
We have no first principles justification for this dependence yet, but the fairness of such 
approximation will be justified later on. 
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Figure 4. The infection probability threshold A* as function of d for typical values of 
the self-isolation probability q of the Q&I process. The lines are fittings based on the 
stretched exponential function A* ~ exp[— (d/a) c \. 



Allowing for the transition values A* obtained from the simulations, we exhibit the 
phase diagram of the model in the plane A versus q for typical values of d in Fig. [5j 
From this plot, it is possible to understand that a disease-free regime is achievable even 
for high infection rates and low compliances with the medical recommendations of home 
isolation, but only when the individuals are near to living alone or the ratio between 
the members of her household by the total number of relationships is small. As the 
fraction of fixed links (relationships) d associated with the household soars, the size of the 
epidemic phase increases as well. This is easily comprehended: since the number of fixed 
links augments, the quenched part of the network becomes more dominant and so the 
number of "channels" through which the epidemics can disseminate. Correspondingly, 
we can end up in a situation for which increasing the consciousness q of the infected 
individual has little or none effect. In this respect, looking at the curve with d = 0.3, 
we verify that the threshold infection rate A* barely varies with the rate q, for q >~ 0.6. 

On the other hand, in the interpretation of this diagram, we can also resort to 
the poll results published by Eastwood and co-workers [TT]. Looking to the data from 
Australia, during the H1N1 outbreak the lowest compliance rate with a public health 
request was of 62.6% for the avoidance of public gatherings. The rate of compliance 
decreases further to as much as 50% when household quarantine is considered [50J. 
As previously mentioned, in the same work of Ref. [50], it was given statistical 
significance to the fact that literacy skills influence the individual performance regarding 
proper compliance with a self-isolation recommendation as well as other public health 
interventions. Bearing in mind that Australia is one of the few countries with a Human 
Development Index (HDI) higher than 0.900, which in total account for less 10% of the 
World's population (World HDI = 0.751), and that it is expected that the lower the 
HDI, the lower the compliance rates, assuming reasonable values for the parameters 
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Figure 5. Phase diagram of the model, separating the disease-free and the epidemic 
phases, for typical values of the density of fixed links d. The symbols represent the 
threshold values A* obtained from the numerical simulations and the lines are plotted 
as a guide to the eye. 

we evidence that it is extremely unlikely that a Q&I policy under such conditions will 
induce a disease- free state. As a matter of fact, for infection rates of approximately 20%, 
as those estimated for the influenza, and d = 0.2 (8 acquaintances per 10 relationships 
and the household composed of three people), we get an epidemic phase if even the 
probability of compliance is 1, which is basically an Utopian value. As regards the 
value of d, we must allow for the fact that at present, people tend to actually (in a 
vis-a-vis sense) interact much less than in the past, particularly in countries where 
internet has had strong implantation. In other words, although the world is apparently 
more connected, activities like home-working and on-line shopping among others have 
boosted their relevancy and with that decreasing the number of acquaintances (and 
casual contacts). 

3. Analytical considerations 

Following the rules of our model, we can do some analytical considerations regarding 
our problem, which corresponds to a network of N individuals and L dual flux links 
that are quenched for the household and allow rewiring for the other cases. At initial 
time, to, the average degree of connectivity is K Q = 2 L / N = (k) and can be split into 

K =K h + K a 

(2) 

= dK + (l-d) K . 
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Assuming all individuals (nodes) equal, we can write a master equation reflecting the 
evolution of the number of infected people which reads, 

I t +i = h + A S t I t - al t 

(3) 

= I t + \(K t -I t ) I t -a I t , 
which in the continuum limit yields, 

^ = X(K t -It) It- al t (4) 

where, 

S t + It = K t . (5) 

When the individual preserves all of its links, K t = Kq (for all t), the stability limit 
of the solution to Eq. (J3J) yielding it_>oo = is given by the infection rate A = a/K as 
given in Ref. [H]. In our case, we must focus on the dynamics of the connections of 
infected people as well. We know that at each time step, an infected person is given 
a medical advice to self-isolate, which implies in a cut of her acquaintances relations. 
In view of the fact that diseased individuals have a probability q of disrupting these 
contacts, the evolution of K a of an infected person is on average given by, 

K? +1 =K?-qK? 

(6) 

the solution thereto is, 

K? = {l-d) K exp[-qt}, (7) 

or, 

K t = K [d +{l-d) exp[-qt]]. (8) 

Now, the determination of the critical threshold from Eqs. (0} and ([6]) can be 
tremendously simplified using the argument [31] that for the epidemic phase to persist 
in time, one must guarantee the infection of other people during her illness span that 
lasts a time of the order of a -1 . In other words, the reproductive rate of infected 
people, R, must be greater than one. At time step t, we can determine the mean-field 
reproductive rate, 

R t = K - [d + (1 - d) exp l-qt]] . (9) 
a 

Averaging R t during the diseased cycle, 

R = a r Rtdt, (10) 
J o 

and solving R = 1 in order to A, we get, 

A c = ~rm 77— ^ ~ (U) 



d M-i -exp -1 (1-d) 
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Figure 6. Left: Critical infection probability of general solution A c vs fraction of 
quenched links d representing the household (symbols) and the numerical adjustment 
considering Eq. (JlJ. Right: Phase diagram of the model considering the general 
solution A c vs the probability q of compliance with a self-isolation recommendation. 
Both curves are obtained from Eq. using the same parameters we introduced 

in the simulations, namely a = 0.2 and Kq = 5. The colour scheme is the same as 
Figs. S] and [5] for easier comparison. 



the limits of which, namely, d — > 1 (quenched network) and d — > (complete 
reconnection with probability q) are verified (T2J [TU [31] . 

In Fig. El we depict A c as a function of d and q, respectively. Comparing both 
panels with Figs. H] and |5l we verify that only the qualitative behaviour of A c concurs 
with that of A*. Specifically, the mean-field solution underestimates the critical values of 
the infection probability. To understand this discrepancy we avail ourselves of recently 
published work [56] where it is ascertained that mean-field like approaches are more 
accurate for networks with high mean degree or high mean first neighbour degree. As 
can be perceived from Eq. flSJ) and using the same parameters values of the simulations, 
the mean degree significantly fluctuates from susceptible to infected nodes. Therefore, 
at the stationary critical state and beyond it, the typical configuration of the network 
includes an important fraction of nodes with a small degree of neighbours, a situation 
that tallies with conditions of inaccuracy we have referred to. 

Nevertheless, we can provide an argument to our empirical adjustment of A* as a 
function of d from Eq. ( II ip . Explicitly, for q > a, that equation can be rewritten as, 

~™ r q] i q 1 1 _1 

Q 



A, 



K (l 

which in first order yields, 
A c ~ exp 



exp 



1 + 



exp 


_q 

a 


a 


- 1 


1 


- exp 


_£ 

a 





d 



d 
V 



(12) 



(13) 
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q 


a' 


b 


x 2 


0.3 


0.44 ±0.02 


0.13 ±0.01 


4.5 x 10~ 4 


0.5 


0.22 ±0.01 


0.25 ±0.02 


4.7 x 10~ 4 


0.8 


0.05 ±0.01 


0.89 ±0.04 


1.7 x 10~ 2 



Table 2. Values of the parameters a' and b to numerically adjust the points of Fig. @] 
with Eq. fll5]). 



with, 



V 



1 - 


- exp 


_£ 

a 




exp 


_q_ 

a 


a 


- 1 



(14) 



A comparison between Eqs. ( TTT1) and ([TBI is presented in the left panel of Fig.[6j Therein, 
for small d the orange dotted line fits the initial points, but then decays faster. The 
stretching exponent c < 1 prolongs the curve. We can also confront our proposal 
Eq. (JTJ) with Eq. ( fill . In Fig. |6] (left side), we have numerically adjusted the points 
given by Eq. (TTT]) with Eq. (Tj[|). As visible, the adjustment matches the curve quite well 
with hardly perceptible deviations. Because of the qualitative similarity between both 
curves, we can be further tempted to consider an adjustment of the data in Fig. H]with 
a mean-field like dependence, 

A * = TT~d]~a'' < 15) 
the values of which are presented in Tab. 2. Comparing the values of \ 2 m both 
tables, we verify that despite the good agreement in the mean field curves, the stretched 
exponential clearly beats Eq. (fT5"]) . especially for large q, which upholds once again our 
previous choice. 



4. Final Remarks 



In this work, we have studied a modified version of the SIS epidemic model that 
takes into account quarantine and isolation practices according to a self-compliance 
with the medical requests. For each person, their relationships are pigeonholed as 
acquaintances and members of the same household. For the former, there is a probability 
that the connection is disrupted, quantifying the degree of consciousness of the individual 
towards the epidemic disease, whereas for the latter the connections are fixed aiming at 
representing the possibility of a compulsory estrangement. For a Susceptible-Infected- 
Susceptible process, our results have shown that quarantines following this structure, 
which agrees to WHO directives, is extremely unlikely to thwart the propagation 
of an epidemic infectious disease using standard figures regarding infection rate and 
compliance with medical (public health) requests. It should be stressed that our results 
must not be interpreted as arguing the uselessness of self-isolation in a overall point of 
view. The outcome of our study indicates that self- isolation is likely to be ineffective 
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per se and signals the importance of coordinated public health or non-pharmaceutical 
interventions in order to control the untamed spread of an epidemic. 

That being said, we would like to highlight that our results can be further explored 
to take into account either different versions considering the impact of other public 
health (non-pharmaceutical interventions) or the joint impact of self-isolation and some 
other intervention(s). Besides scrutinising different types of epidemics, we can also study 
the problem of isolation using mixed topologies [57] for the social network which depend 
on the type of relationship between people and institutions. Regarding this point, the 
problem can be honed by bolstering the acquaintance type connections, i.e., in the 
rewiring step people are biased to rewire to acquaintances and members of the same 
household of people with whom they already have a relationship. Another important 
investigation concerning the dynamics respects the assumption of an (average) isolation 
time instead of a fixed (unitary) time scale we have taken into consideration. 
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